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Abstract 



^^. This paper is dedicated to the study of the existence and the properties of electron- 

(— I i electron bound states in QED3. A detailed analysis of the infrared structure of the 

perturbative series of the theory is presented. We start by analyzing the two-point 
Green's function, in the Bloch-Nordsieck approximation. The theory appears to be 
plagued by severe infrared divergencies, which nevertheless disappear when vacuum- 
ed I polarization effects are non-perturbatively taken into account. The dynamical induc- 
tion of a Chern-Simons term is at the root of this mechanism. 

From the inspection of the electron-electron non-relativistic potential it then fol- 
lows that equally charged fermions may either repel or attract and, moreover, that 
bound states do in fact exist in the theory. We calculate numerically the binding ener- 
gies and average radius of the bound states. We find an accidental quasi-degeneracy of 
the ground state of the system, between the lowest-energy I = —3 and / = —5 states, 
which could be related to a radio-frequency resonance in high- Tc superconductors. 



* Supported in part by Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq) and by 
Fundagao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP), Brazil. 

^On leave of absence from Instituto de Fisica, Universidade Federal do Rio Grande do Sul, Caixa Postal 
15051, 91500 - Porto Alegre, RS, Brazil. 

■fOn leave of absence from Instituto de Fisica, Universidade de Sao Paulo, Caixa Postal 20516, 01498 - 
Sao Paulo, SP, Brazil. 



1 Introduction 

As is well-known [|[], the Maxwell-Chern-Simons (MCS) theory is a (2 + 1) -dimensional 
model describing the coupling of charged fermions {ip, ip) of mass m and electric charge e 
to the electromagnetic potential A^ via the Lagrangian density 

z — % — — — 

+ -tljYd^,'ip - -{d^^)Y4j - m^^ + eip-i^'A^ip, (1) 

where F^y = d^Ay — dyA^, 9 is a. topological parameter with dimension of mass and A is 
a gauge-fixing parameter. Neither parity nor time reversal are, separately, symmetries of 
the modelQ. The lowest order perturbative contribution to the effective fermion-fermion 
low-energy potential, arising from (0), was recently computed and reads 



V(r) = — 1 \Ko (\e\r) - — [1 - \e\rR\ (\e\r)] L. (2) 

Here, r is the relative distance between the electrons, r = \f\, L = rip2 —r2Pi is the orbital 
angular momentum, whose eigenvalues are denoted by /, p is the relative linear momentum 
of the electrons, and Kq and Ki designate the modified Bessel functions^]. 

The linear dependence of V^ on L accounts for the breaking of parity and time reversal 
invariance in the non-relativistic approximation. The aim in Ref. p| was to determine 
whether the potential (^ could bind a pair of identical fermions. For positive values of 6, 
a numerical solution of the Schrodinger equation confirmed the existence of a bound state 
for e^ /{ttO) = 500, m/6 = 10^ and 1 = 1. Further numerical analysis indicated that all 
identical-fermion bound states are located in the region e^ /{t^O) > 1. 

Thus, we were naturally led to consider the limit 6' — > where the MCS theory degen- 
erates into QED3, whose Lagrangian is 

C = --J^yF^^- - lj^d^A^){dyA'') + '-^i^d,^ - '-{d,^)ri' - m^^ + e^YA,^- (3) 



^Throughout this paper we use natural units {c = h ^ 1). Our metric is goo = ^511 — ~522 = 1, while 
for the 7-inatrices we adopt the representation 7^^ = (73, 7^ = ia^ , 7^ = ia'^; ct', i = 1, 2, 3 are the Pauli spin 
matrices. 

^The potential (||) looks similar to the one derived by Kogan [Q for the same problem. There are some 
important differences, however, that we would like to stress. The derivation in Ref. is solely based 
on relativistic quantum field theory Q. On the other hand, in Ref. |^ the corresponding potential is 
determined in two steps. First, the A^ vector created by a point charge is computed. Then, the quantum 
mechanical Hamiltonian describing the low energy relative motion of the two electrons is assumed to be 
that of a charged particle in the presence of the external field A^. As a consequence, an extra term 
proportional to [1 — |6'|riiri (l^lr)] arises in the formulation of Ref. 0. 



Power counting indicates that the S-matrix elements of QED3 are plagued with infrared 
singularities whose degree of divergence grows with the order of the perturbative expansion. 
The situation is, therefore, more severe than in QED4 where all (on-shell) infrared diver- 
gences are logarithmic. Hence, the usual mechanism (soft bremsstrahlung) for cancelling 
infrared divergences at the level of the cross section fails, since it would only be opera- 
tive for the leading divergence. Furthermore, a non-perturbative analysis of the two-point 
fermion Green function G{p), carried out within the Bloch-Nordsieck (BN) [^-|^ approxi- 
mation, shows that G is well defined for generic values of the momentum p but develops an 
essential on-shell singularity. Through a multiplicative renormalization, G -^ Gr = Z^^G, 
one may define a new Green function behaving on shell as a simple pole. However, Gr is 
found to exhibit an essential off-shell singularity which reintroduces the difficulty for the 
scattering amplitudes. All these developments are presented in Sec. 0. 

We recall next that the BN approximation does not take into account vacuum polariza- 
tion effects. This does not constitute a serious limitation for QED4, as far as the infrared 
structure of the theory is concerned. In fact, for QED4 these diagrams are, on the one 
hand, free from infrared singularities while, on the other hand, they do not give rise to 
a mechanism for dynamical mass generation. In contrast to this, when vacuum polariza- 
tion effects are non-perturbatively incorporated into the photon propagator, the infrared 
behavior of QED3 changes drastically, as is demonstrated in Sec. ^. Essentially, the the- 
ory is reformulated in terms of a massive vector boson whose mass is \6\ = e^/(87r), the 
dynamically induced Chern-Simons term |1|,^ being at the root of this mechanism. Thus, 
the infrared disease affiicting perturbative QED3 is cured and one is allowed to read off the 
corresponding effective electron-electron low-energy potential directly from Eq. (^, after 
replacing 9 — > — e^/(87r). Hence, 



eV^QED3(^-^ 
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(4) 



The terms proportional to Kq in Eq. (^ are both repulsive. The term SL/lmr"^) becomes 
attractive (repulsive) for negative (positive) eigenvalues of L, while the term proportional 
to Ki acts in the opposite way. Hence, unlike the case of QED4, electrons may repel or 
attract in QED3, the infrared structure of the theory being responsible for this unusual 
behavior of the interaction between equally charged particles. A numerical investigation 
on whether the potential (^) can sustain electron-electron bound states is presented in 
Sec. |[ There we show that such bound states do in fact exist. For some of them the 
corresponding dissociation temperature will be seen to be in qualitative agreement with 
the values encountered for the critical temperature of high-Tc superconductors. 

Our conclusions are contained in Sec. ^. A preliminary version of some of the results 
presented in this paper is contained in Ref. |T^ . 

We remind the reader that the signs of the electron's spin and of the induced topological 
mass are both equal to the sign of the fermionic mass term in the Lagrangian (^) |I|]. Thus, 
in a theory involving both spin-up and spin-down fermions there will be no induced Chern- 



Simons term and, for the specific case of electrodynamics, there will be no cure for the 
infrared problem. This explains why it is essential that we restrict ourselves, in this paper, 
to the case of a single two-component fermion field. The generalization to the case of 
several flavors is straightforward and will be mentioned along the way. 

2 Infrared Structure 

As we already mentioned in the Introduction, power counting alone reveals that perturba- 
tive QED3 is afflicted by severe on-shell infrared singularities. To gain some insight into the 
infrared structure of this theory, we shall start by computing the two-point fermionic Green 
function in the BN approximation. To keep the singularities under control, we approach 
QED3 as the 6-^0 limit of the MCS theory (|l]). Clearly, both theories exist off shell and 
MCS goes into QED3 as ^ 0. 

The BN approximation consists of the replacement of the 7'^-matrices in Eq. (|1]) by a 
vector u^, with v? = 1. As a consequence, the usual free fermion propagator is replaced by 
the retarded function 
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where H(x^) is the Heaviside step function. Then, the interaction does not imply any cor- 
rections to the vector meson propagator. For our present purposes this does not represent 
a serious drawback since vacuum polarization insertions do not alter the leading infrared 
behavior of a graph. 

The Green functions of the theory will be computed by functionally differentiating the 
generating functional Uo[J^, r], fj] with respect to the external sources. In this paper, J^ and 
?7, fj denote the vector meson and fermion external sources, respectively. After integration 
on the fermionic degrees of freedom one finds 

+ i I d^x Jf,{x)A^{x) -i I d^x I d^y f]{x)G[A\x, y]r]{y) \ , (6) 

where e^ is the fermionic determinant, S[A] = —F^^F^'^ /A and G[A|x, y] is the fermionic 
two-point Green function in the presence of an external field A^. 

The solubility of the model in the BN approximation is partly due to the fact that 
the A^ propagator is not corrected by the interaction, and partly due to the factorization 
property of (^[Alx, y\. In fact, one can convince oneself that the differential equation 



is solved by 



where 
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G[A\x, y] = h[A\x]GF{x - y)h-'[A\y], 



h[A\x] = exp 
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One can check that expression (^ can be cast in a form analogous to that quoted in Ref. [§ 
for the case of QED4. 

Formally, the computation of D yields D = J^ de' J d^xG[A\x, x]u^An{x), which in view 
of d^) reduces to 

D = e d^xGF{x,x)u^Af,{x), 

showing that D can, at most, depend linearly on A^. However, as seen from (j^), Gf{x.,x) 
is ambiguous. In the original model, with 7^^ instead of u^, Lorentz invariance demands 
the vanishing of the tadpole contribution to the fermionic determinant. Based on this fact 
we shall therefore take, from now on, D = 0. 

The complete two-point fermion Green function can now be readily found. 



G{x,y]m,9) = 



1 lS^Uo[J, = Q,v,v] 
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where 

S^{x,y-z) = i^^{x-z)~i>^{y-z), 
and il^^^ is the photon propagator. From (|^) it follows that, in momentum space. 
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Hence, for the Fourier transform G{p; m, 9) = j d?xeiq){ipx)G{x\ m, 9) one gets, after some 
algebra, 

/•oo 

G(p; m,9) = / duexp \iv{u ■ p — ni + it) + /(z/, 9)] , 

JO 

where 



and 



/(^, 0) = -j^, I d'k{^ - cos[z/(n ■ k)]}U{k, 9), (9) 



U{k,9) 



{u ■ kf 



Since u^e^^pU'^ = 0, the second term in (|) does not contribute to U{k, 9). Thus, the 6' — > 
hmit of G in the MCS theory is the same as that of the corresponding function in massive 
QED3. 

The behavior of the integrand in Eq. (|^) in the hmit k —* 00, with u ■ k kept constant, 
tells us that f{i>, 9) develops a logarithmic ultraviolet divergence. This divergence, which 
can be absorbed into a mass renormalization, is entirely due to the approximation being 
used {u^ instead of 7^)0. The two-point Green function G in terms of the renormalized 
mass mji turns out, then, to be given by 



G{p; niR, 

where Q = u ■ p — rriR, Ei{z) is the exponential-integral function, 

^n = m+-ln^-j, (11) 

A is an ultraviolet cutoff, and M plays the role of a subtraction point. For 6' 7^ 0, the 
singularity structure of G{p; uir, 9) can be seen to consist of a pole alQ 

u.p = mR + — [ln—-\ 
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■^A similar situation arises in QED4 where the corresponding function f{v) presents linear and logarith- 
mic ultraviolet divergences, the linear divergence being induced by the replacement 7'^ ^ u^ [^. 

"^An analogous situation arises in the one-loop calculation of Ref. [|l|, where a gauge-dependent contri- 
bution to the electron mass is found. 



plus branch cuts. If one insists in keeping the pole of G at the renormalized mass defined 
in Eq. (|ll|), one must choose M = 9exp{\), thus lumping all the gauge dependence of G 
into the subtraction point. Moreover, the Green function G can be renormalized, 

G{p; mn, 6) -^ Gr{p; tjir, 9) = Z~^G{p] mn, 9), 
such that 

lim {u-p-mji)GR{p]mR,9) =i. 

It may be checked that the ultraviolet-finite renormalization constant 

does the job. Everything runs smoothly for 9 y^ 0. 

However, the above scheme does not survive the limit 9^0. While the pole of Gr 
can be kept at the value of the renormalized mass, the constant Z presents an essential 
singularity which in turn is transferred to Gr and shows up for all off-shell values {u ■ p ^ 
niR) of the momentum. One may try to circumvent this difficulty by defining the two-point 
fermion Green function of QED3 as the 9^0 limit of (p!0D , keeping M arbitrary. For Q < 
and after a Wick-like rotation one arrives at 



I r°° , I e^ 
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where C is the Euler constant. On the other hand, if Q > one obtains 
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although in this last case the Wick-like rotation is only allowed if e^ 7^ 0. Expressions (|I^ 



and (|T3D are well defined for generic values of Q but contain an essential singularity at 

In both cases above, the BN approximation fails to provide an acceptable functional 
form for the two-point fermion Green function of QED3. This is an indication that the in- 
frared singularities occurring in each order of perturbation theory do not add up to a finite 
limit. Moreover, the inclusion of soft bremsstrahlung emission, which in QED4 removes the 
infrared problem from the cross sections, is of no help in the present case since it would 
only act on the leading divergences. Nothing changes, qualitatively, if vacuum polarization 
insertions are added order by order in perturbation theory. Nevertheless, when vacuum po- 
larization effects are non-perturbatively incorporated into the photon propagator, a photon 
mass is dynamically generated and as a consequence the infrared disease is cured. To show 
how this comes about is the purpose of the next section. 

^It is instructive to compare this result with the corresponding one in QED4, where the two-point 
fermion Green function exhibits a power law behavior in the variable Q ra. 



3 Vacuum Polarization Effects 

The lowest-order graph contributing to the vacuum polarization tensor 11''°" yields 

^'^' '^ J {2tt)^ [(g + fc)2-m2](P-m2) ' 

Since we are interested in the quantum corrections to a non-relativistic potential, we shall 
retain only those terms which are zero and first order in q. Gauge invariance alone ensures 
that n^°'(0) = 0. For the first-order contribution, which gives origin to the induced Chern- 
Simons term |]I],3, one finds 

We emphasize that 11^^^ (g) is ultraviolet-finite and that, therefore, no regularization is 
needed for its computation. Alternatively, if one adopts the point of view that ultraviolet 
divergences should be kept under control by the regularization of the entire theory, the result 
quoted in this last equation is true only if a parity-time reversal invariant regularization is 
used. 

We now modify the photon propagator by resumming the geometric series resulting from 
the iteration of li^pj (see Fig. [l|). This is formally equivalent to adding to the Lagrangian 
(H) the term 

which is nothing but a Chern-Simons term with and induced topological mass 

0^n = -f- (14) 

Hence, in the low-energy regime, the theory effectively describes the interaction of fermions 
with massive vector particles of mass \6in\ and is, therefore, free from infrared problems. 
Moreover, from Eqs. (0) and (|1^ it follows that the QED3 effective electron-electron low- 
energy potential is, in fact, the one given in Eq. (^. 

We verify the consistency of the present formulation of QED3 by showing that all 
remaining contributions to V^^^^ are, up to some power of \n{e'^/m), of order e^/m or 
higher with respect to (^, and vanish as e^/m — ;► 0. To see how this comes about, we 
compute the vertex correction A^ to V^^^^ arising from the diagram in which only one 
massive vector particle is exchanged^ (see Fig. |^). After the replacement 6 -^ 6in, the 
massive vector field propagator can be read off directly from Eq. (4) of Ref . , 

^For the vertex corrections in MCS theory see Ref. [0. We would hke to point out, however, that the 
hmit — > in the MCS theory differs radicaUy from the hmit 6^ = — e^/(87r) ^ in QED3, since in the 
latter case the vanishing of the topological mass also implies the vanishing of the coupling constant. 
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where P^v = g^lv — k^k^j/k"^ , and an arbitrary function /(fc^) has been incorporated into the 
longitudinal part. Accordingly, A'^ can be split as follows 
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where the subscripts g, e, and L make reference to those pieces of D^^ proportional to (7^,^, 
^fiupk'^ and kf^k^, respectively. The computation of A^, is straightforward and yields 



A1 = tXv''+\p[)Yv^-h 
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where, as in Ref. 0, f '•^-'(pi) (■^'^'^•'(p']^)) is a two-component spinor describing a free elec- 
tron of two-momentum pi (p'^) in the initial (final) state. Since A^ does not depend on the 
momentum transfer q = p'l—pi, it can be absorbed into the wave-function renormalization 
constant. When computing A^ and A^ only zero- and first-order terms in q/m will be 
retained, since our interest here is restricted to the non-relativistic regime. Furthermore, 
all the momentum integrals over loops are ultraviolet-finite and there is no need for reg- 
ularization. After absorbing the zero-order terms into the wave-function renormalization 



constant, one arrives at 
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From Eqs. (IHHIl) it follows that A^*-^''(g) behaves as (e^/m^jln^e^ /m) when e^/m -^ 
{6in/m -^ 0), while power-counting indicates that the full vertex insertion A'^(g) (Fig. |^) 
may diverge logarithmically at the just mentioned infrared limit. The fact that A'^''^^(g) 
exhibits an improved infrared behavior is not a peculiarity of the particular insertion under 
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analysis but applies to any vertex part involving an arbitrary number of exchanged massive 
vector particles. Indeed, the leading infrared divergence of any of these parts only shows 
up in those terms containing even powers of the momentum transfer g, as can be seen 
by setting to zero the loop momenta in the numerators of the corresponding Feynman 
integrals. Thus, the terms linear in q exhibit a milder infrared behavior. 

By using the technique described in Ref. one finds the correction Al^^^^^ arising 
from the diagrams in Fig. ^, 
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From Eqs. (^^-|T^ it follows that eAy*^^°^ also behaves as (e^/m) ln(e^/m) when e^/m — >■ 
0, and turns out to be negligible if e^/m <^ 1. This establishes the region of validity of our 
approximat ions . 

4 Electron-Electron Bound States 

We turn next to the investigation of the existence and properties of bound states of two 
identical fermions of mass m and charge e interacting through the non-relativistic potential 
yQEDs gjygj^ jj^ ^g) The corresponding radial Schrodinger equation is found to read 

HiRniiy) = GnlRnliy), 

where Rni is the radial part of the wave function, / is, as defined before, the angular 
momentum eigenvalue, n is the principal quantum number, and Tii is the effective radial 
Hamiltonian, namely. 



HiRniiy) 



d^Rniiy) ^IdRniiy) 



+ Uf^iy)RUy)- 



dy"^ y dy 

Here t/f denotes the effective radial potential arising from (^, 

Uf\y) = 4 + ^(1 + /5-) Ko{y) + ^ [1 - yK,{y)] , (20) 

y 2 y^ 

where we introduced the dimensionless parameters y = e^r/(87r), ain = e^/(7r|6'j.„|) = 8, 
Pin = rn/\6in\ = STcm/e'^, and e„; = Gin'^mEni/e'^, Eni being the energy eigenvalue. 

A straightforward analysis of (pOf ) reveals that electron-electron bound states are pos- 
sible only for / = — 1, —3, —5 and —7. In order to see this, we consider the derivative 
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For I negative and outside the interval [—7,-1], this derivative is negative definite and, 
therefore, t/f has no local minima. Inside this interval, on the other hand, an examination 
of the behavior of f/f reveals that for small y it diverges to +00, while for large y it 
approaches zero from negative values. It then follows that t/f has an absolute minimum, 
thus opening up the possibility for bound states. The fact that Fermi statistics requires 
that / be odd reduces the possibilities to the values mentioned above. Had we chosen the 
negative sign for m we would have obtained 9in = +e^/(87r) and bound states, with identical 
energy eigenvalues, would be possible for / = +1, +3, +5 and +7. This analysis can be 
easily extended to the case of two different flavors. If the electrons are in an antisymmetric 
flavor state, bound states may emerge for / = —2, —4 and —6. 

The straight numerical resolution of the Schrodinger equation is useful for determining 
quickly whether bound states exist or not, but it is not so useful for the determination of the 
binding energies of such bound states. The reason is that in order to solve the equation one 
needs, beforehand, a deflnite value for the energy parameter. Since the energy eigenvalue 
is not known a priori, one must search for it, solving the equation repeatedly by starting 
with the boundary condition for Rni at y = 0, and propagating the wave function forward 
until one realizes the boundary condition for large y. Besides being a difficult process in 
itself, due to the very sensitive dependence of the asymptotic behavior of the wave function 
on the value of the energy parameter, it interacts in a complex way with the finiteness of 
the integration interval and of the integration range. 

The existence of bound states for the potential given in (^) was examined numerically 
by means of a stochastic variational algorithm. We found this method much more con- 
venient than the resolution of the complete eigenvalue problem, which is large and hard. 
The variational algorithm consists of a direct search for the lowest eigenvalues only, and 
is therefore more appropriate for our purposes. In all cases which were examined with 
both this procedure and the Schrodinger equation, the results were the same, within the 
numerical limitations of each method. With the variational algorithm we were able to 
quickly identify, for a given /, the state of minimum energy (n = 0), for a variety of values 
of Pin, and to estimate roughly the binding energies. However, for a precise determination 
of the binding energies, long double-precision runs were necessary. We also measured the 
expectation value (f„/) of the radius in the state thus obtained. A typical configuration of 
Rqi, and the corresponding potential, are displayed in Fig. |[ 

In the variational algorithm we varied the wave function Rqi, subject to the boundary 
condition Roi{0) = 0, while trying to decrease the energy of the state, which is given by 
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During the variation of the wave function, besides implementing the boundary condition 
at the origin, we kept constant its maximum value. This implies that the normalization 
of the wave function relaxes to some arbitrary value. The requirement that the wave 
function should decrease exponentially for large values of y was realized automatically by 
the algorithm. 

In order to be able to represent the wave function R^i and the expectation value of the 
energy e^i numerically, we must chose a discrete collection of points (sites) along the radial 
direction. Hence, we calculated the ratio of sums 
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El=r Ay y^ Rl, 



where ARqi are forward differences, and the values of Ay and imax were chosen to ensure a 
big enough range of integration Y = Ay imax and a small enough integration interval Ay. 

We then varied the wave function randomly, accepting any changes that decreased the 
energy. The function was varied by sweeping the sites, excluding the point at the origin 
but including the point corresponding to F, where it was left free to fluctuate. At each site, 
the wave function was changed by adding to it a random number in the interval [— £,e). 
After a certain number of sweeps (blocks), the position of the maximum was determined 
and the function was renormalized to keep constant its maximum value. 

During the sweeping procedure we monitored the acceptance rate of the trial changes in 
the wave function, which depends strongly on the value chosen for the range e. In order to 
keep the code efficient, it is necessary to ensure that the wave function is being changed as 
fast as possible. The optimum value for e is the one which maximizes the product of this 
range by the acceptance rate. This product is a good estimator for the overall speed with 
which the wave function is being modified. We included in the code a feedback mechanism 
which searched for this optimum value of e by maximizing the product. 

The runs were started with a large value for e, which was then continuously adjusted 
by the code as the run proceeded. We found that the feedback mechanism worked very 
well to decrease e to appropriate values but that sometimes it could overshoot and get 
hung at extremely small values of the range. This was particularly true in the long runs 
needed for the precise determination of the binding energies. This problem was overcome 
by a procedure which we call quenching. We stopped the runs two or three times, storing 
the wave function produced so far, and restarted them with a large value for the range e. 
This procedure greatly increased the convergence speed of the code. Figs. ^ and |^ display 
the relaxation of the binding energy and of the average radius of the bound state, at the 
final quenching cycle, in a particular case. 

Another relevant numerical issue is a phenomenon related to the finite numerical preci- 
sion of the computer. During the sweeping procedure it was often found that many changes 
in the wave function did not produce any changes at all in the energy. Clearly, in these 
cases the change in the energy was below the numerical precision available. When this 
happened we accepted the change anyway, monitoring separately the acceptance rates due 
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to relaxation of the energy and to this phenomenon. The acceptance of these changes in- 
teracted with the relaxation changes, increasing the efficiency of the code. Another source 
of numerical errors was related to the sweeping procedure. In order to save computation 
time, during the blocks of sweeps, the energy was not re-calculated every time a site was 
upgraded, but had instead its variations added to it along the way. This introduced numer- 
ical errors in the system and, in order to keep the errors small, after each block of sweeps 
the energy was calculated anew. 

The finite-precision phenomena described above caused a small random drift in the 
results, which were used to estimate the stochastic errors. The stability of the results 
was tested against variations of Y and Ay. Due to the exponential decay of the wave 
function for large y, it was found that, given a large enough Y, the results were essentially 
independent of it. They were somewhat more sensitive to the finiteness of Ay, approaching 
a limit linearly with Ay, as Ay -^ 0. Also, the code became progressively less efficient as 
one decreased Ay. The reason for this is easily understood. For small Ay any change in 
the wave function at a single site tends to produce large spikes in its derivative and, then, 
the derivative term of the energy causes all but very small changes to be rejected. For 
small Ay the local upgrading is inefficient to change the long-wavelength components of 
the configurations. 

Being prevented by the computational limitations from decreasing Ay indefinitely, we 
estimated the Ay = limits of the results by linear regression to zero from three finite- Ay 
runs, for each set of values of the parameters j3in and /. These are the data presented in 
Figs. [7|-[TI|. They were calculated assuming m to be the usual electron mass. The energy 
eigenvalues for / = —5 and / = —7 are very close to those corresponding to / = — 3 and / = 
— 1, respectively, and for this reason they have not been included in the figures. The error 
bars are too small to be visible in the graphs. The stars mark the points calculated directly, 
while the lines were obtained by splines interpolations. The corresponding dissociation 
temperatures Ta^ni are plotted in Fig. ^ Table |l] contains approximate fits for the data. 

An examination of Ui in Eq. (^) shows that for large y the potential changes very 
little between the cases / = —3 and / = —5, and between the cases / = —1 and / = —7. The 
dominant term /(/ + 8)/y^ does not change at all; the only term that does change contains 
Ki, which is exponentially damped for large y. It is not so clear how much the potential 
changes near the well, which is the most important region as far as the bound states are 
concerned. In fact, one finds but very small differences between the energies for / = —3 
and / = —5, and between those for / = —1 and / = —7. The ground state of the system 
is the / = —3, n = state, but the / = —5, n = state is very close to it, forming a 
quasi-degenerate doublet. 

Assuming that transitions between these states do occur, the transition energy would 
correspond to very low resonant frequencies, in the radio spectrum. This resonant frequency 



/ij appears in Fig. |T0|, as a function of Pin, and in Fig. |TT|, as a function of the dissociation 
temperature. Presumably, the states with larger n also form similar doublets and there 
should be further transition frequencies, similar to this one, corresponding to them. 
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5 Conclusions 

There are ranges of the parameters / and (3in where the results of Sec. ^ are numerically 
consistent with the observed phase-transition temperatures of high-Tc superconductors. 
We do not think that this fact is merely accidental, although we are not claiming that 
QED3 appropriately describes all features of high-Tc superconductivity. Also note that 
if the initial state of the system contains only low energy electrons, energy conservation 
forbids the production of electron-positron pairs. No electron-positron bound states can 
be formed and, as a consequence, only the electron-electron bound state can characterize 
the many-body ground state. 

Certainly, the calculations presented here are still incomplete. For instance, one should 
recalculate the potential for the finite-temperature field theory, which presumably will 
introduce corrections in the energies and temperatures. It also remains to be shown that 
the theory exhibits a phase transition and that the many-body ground state is a condensate 
of electron-electron pairs. However, if this turns out to be the case, we find reasonable to 
think that the dissociation temperatures found in this paper should be, at least, rough 
approximations for the true phase-transition temperatures. 

We believe that the accidental quasi- degeneracy uncovered in this work may be impor- 
tant to test the relevance of the model for the description of high-Tc superconductors. The 
corresponding transition frequency is directly related to the binding energy of the ground 
state, and therefore to the predicted dissociation temperature. The relation thus implied, 
between this frequency and the phase-transition temperature, may provide a means of test- 
ing the applicability of the model for the description of high-Tc superconductors. If the 
model considered here is more than superficially related to high-Tc superconducting mate- 
rials, it should be possible to detect such a transition. It certainly would be very interesting 
to further pursue this matter. 
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FIGURES 

1. Wavy lines represent free photons, while dashed lines refer to massive vector particles. 

2. The vertex insertion A^. 

3. Diagrams contributing to the potential AV^^^'^. 

4. The bound-state wave function Rqi (dashed line) and the dimensionless potential f/f 
(continuous line) for / = —3, /9j„ = 2000, Ay = 0.25 and Y = 200. The vertical scale 
refers only to the potential. The normalization of the wave function is arbitrary. 

5. Relaxation of the dimensionless energy parameter eo; for / = —3, /9j„ = 2000, Ay = 
0.25 and Y = 200, on the last quenching run, as a function of the blocks of sweeps. 
The vertical scale is relative to the minimum value eo-3 = —0.03266121075526209, 
which corresponds to the zero of the scale. 

6. Relaxation of the dimensionless average radius yoi for / = —3, Pm = 2000, Ay = 0.25 
and Y = 200, on the last quenching run, as a function of the blocks of sweeps. 
The vertical scale is relative to the mean value yos = 18.02868690614052, which 
corresponds to the zero of the scale. 

7. Binding energy Eqi, in eV, versus the dimensionless parameter /9j„, for / = — 1 and 
/ = -3. 

8. Average radius f^i, in A, versus the dimensionless parameter /9j„, for I = —1 and 
/ = -3. 

9. Dissociation temperature Td^i, in K, versus the dimensionless parameter (3in, for 
/ = —1 and / = —3. 

10. Resonant frequency J'r, in MHz, versus the dimensionless parameter /9j„, for the 
(/ = -3)— (/ = -5), n = doublet. 

11. Resonant frequency /r, in MHz, versus the dissociation temperature T^,oz, in K, for 
/ = -3. 
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TABLES 



1. Approximate fits for tlie numerical data. 



Grapli 


Fits 


Eqi X Pin 


Eo_i ^ 2.25 X 10^ A-2-224 eV 


Eo-3 ^ 9.55 X 10^ /5-2-229 eV^ 


foi X (3in 


fo-i ^ 4.05 X 10-2 /^i;ii3 I 


fo_3 ~ 2.86 X 10-2 /3,?„ii7 ^ 


Td,Ol X Pin 


T,,o-i ^ 2.61 X 10^ A-2-22^ K 


T,,o_3^1.11xl09An'-'''^ 


Ir X Pin 


fn ^ 1.29 X 10^3 /5r3-465 mHz 


Ir >< Tdfl-3 


fn ^ 1.13 X 10-1 nf4{K) MHz 



17 



